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ABSTRACT 

Analyses show that the central place in the mathematical description of the spatial-temporal dynamics of one or 
several interacting populations take the task of research of interaction of the migration population processes with 
demographic. Population models that take into account only the demographic processes are well known and fairly 
well-developed [3,4]. This so-called point model, the main assumption which is the assumption of "infinitely fast" stirring 
of individuals in a given area. This assumption is true if modeled area is quite small compared with the mean free path of 
individuals, or, equivalently, with a radius of individual activity [3]. If this provision violated, in population models must 
account migration. 

KEYWORDS: Biological Population, Nonlinear Model, Parabolic Type Equation, Self, Similar Solutions, Lower and 
Upper Solutions 

INTRODUCTION 

The simplest and most widely used in the present situation is the hypothesis of randomness "walk" of individuals 
in space. This assumption allows substantiating of using as a tool of modeling equations of reaction-diffusion type, where 
as a reaction part was used the right part of the point models, and diffusion coefficients (mobility of individuals) are 
assumed to be constant. 

In the framework of these models it is possible to explain such effects as wave propagation of population at 
settling the area and existence of complex spatio-temporal dynamics of population. 

Consider in Q T = [0, T]xR N generalized reaction - diffusion task of Kolmogorov - Fisher type in the 
following form 
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which describes the process of biological populations with double nonlinearity, which diffusion coefficient is 

p-2 

. Here D, m, p, /B > 1 - are given constants, 
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Equation (1) is a generalization of the simplest diffusion model for the logistic model of population growth [1]. 
It can be regarded as the equation of nonlinear filtration, thermal conductivity in the presence of simultaneous influence of 
a source and acquisitions. 

This equation in the case of m=0, k(t)= k >0 constant qas suggested (as well as the cubic equation instead 
quadratic nonlinearity in the right side) Fisher [3] as a stochastic model of the deterministic version of the gene in a 
favorable distribution of diploid populations. He examined in detail the equation and obtained some useful results. 
Heuristic and based on genetics conclusion of the equation have also led by A. N. Kolmogorov, I. G. Petrovskii and N.S. 
Piskunov, classical work [4] which served as the basis for a more rigorous analytical approach to the Fisher equation. 

Properties of the solution of the Cauchy problem and the solution of equation (1) in the case of a homogeneous 
environment, when (1) m = 0 and k(t)=k - constant studied in detail by many other authors (see for example [3], 
where you can find links to other work). Suppose that in (1) m = 0, k{t) = k -constant. In this case, for speed of 
propagation of the wave with wave solutions of (1) 



authors of [4] were rated an estimate C > 2 V kD . 

With regard to the properties of solutions of the initial task for the equation (1) is unknown and so it is interesting 
to trace the evolution of the process of reaction diffusion in heterogeneous environment and explore the impact of 
heterogeneity and case of dependences of the reaction coefficient from the time, i.e. k= k(t). 

This paper is devoted to the study of properties of solutions of the task (1), (2). Clarified conditions on parameter 

m and diffusion coefficient k (?) , where take the place finite speed of movement of a wave. Proved also, bilateral 
estimation of a solution of the task (1). 

We construct self-similar equations for (l)-(2)-more simple for research equations. We construct self-similar 
equation by the method of nonlinear splitting [1] 

Below we propose an algorithm based on the splitting of the original equation of parabolic type with which 
assesses estimations of lower and upper solutions of the task (l)-(2). 



Df'-cf'+f(\-f p ) = Q 



u 



{t, *)=/(£), % = ct + x, 




ALGORITHM OF SPLITTING OF THE PARABOLIC TYPE EQUATION 



We solve the equation 



U t = ku (1 — W ) , 



(3) 



and then search solutions of the equation (1) in the form 
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U(t,x) — U (t)w(T,x) 

Then supplying (4) into (1) we have again equation of the form (1) 

p-2 

m 1 



(4) 
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+ k l (t)w(l-w P ). 



(5) 



where k x (t) = k(t)u" ,(p - 2)+m = k(t)u' n[p - l) . 
only difference is that instead of the t stands T(t) = j" 



J3e 



m(p-2)kt 



m(p-2) 



(l + e mip ~ 2>kl ) 



dt , and instead of k(t) - k i (t) . 



Hence, if mp — m = 0 equation (5) has exactly the form (1), so again, we get an equation of the form (1) 
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dt dt 
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(6) 



m(p-2) 



however, instead of a variable by the time t stand the function t(t) = j* \u (t)] P dt . 
Now consider equation 

OW — m „ i — m 

— = V(IVw \"~ 2 Vw ) 
dr 



(7) 



Equation (7), which is called the etalon for (1), has six types of self-similar solutions, one of which has the form: 

*) = /(£), Z=\x\/T Vp (8) 
where f(£) satisfies the equation 
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(9) 



at £ < (a / b) ( p 1)1 p . From (6) it is clear that if selected as the comparable function z(t, x) — u(t)w(T, x) and 
W<1 in Q T , where W is solution of (6), lt(t) given by (3), it will be a lower solution of the original task, if 
U 0 (x) > z(0, x) . This result can be improved by applying again algorithm of nonlinear splitting applied above. 

Now after putting w(T, x) = f(<^), ^ =\ x \ I T 1 p , equation (6) will turn into a self-similar equation 
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Now if we construct a function z(t,x) — u(t)w(t, x) , where w(f,x) - upper solution of the etalon equation 
(7), for which runs conditions Az < 0 in Q T and U Q (x) < z(0, x) , then function z(t,x) will be an upper solutions of 
the task(l), (2). Here 



Au 



du d 
dt dt 
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du' 
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' ^ IP 

du 
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+ k(t)u(\- U P). 



On these considerations are based algorithm for the construction of upper and lower solutions of the Cauchy task 

(1), (2), the global solvability of the initial task and definition of the critical value of the parameter P , where solution of 
the task becomes unlimited, and changes view of the asymptotic representation. 

T , ■ Q =Q, £>={(t,x): t>0,lxl</(0} i /<7)>0 t>0 

Let's introduce notations ^' lv ' WJ , where lyl) ^ u at ' ^ u - continuous 



function. 



Takes place the following lemma of comparison of solutions [3,5]. 



Lemma. Let X ^ > ^ - generalized solution of the task (1), (2) and functions u ±(^ x ) sucn that u ± > ^ j r 



u =0 u + eC'' 2 (D)nC(D) 

M ± u in Q\D, ± ' ' x v ' ' . 



du' 



dx 



p-2 



dx 
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and in D runs an inequality 



Au + < 0, Au_ > 0, u 0 (x) < w + (0, x), u 0 (x) > w_(0, x) . d n u 0 (x) < u + {t, x), u 0 (x)>u_(t,x), 

111 . 1 I lL 1 1 

inQ. 

Functions U + , U_ are called, respectively, the upper and lower solution of (1), (2). 
By direct calculation we can ensure that the function 

_ p-i 

( — Y 

/„(£)= C\-h\gV 1 (ID 



where Cj < 0, b — 



1 m(p -1) -1 

a , is generalized solution of the equation 



m 



, e df 



+ a^^- + aNf = 0, 



where constant a > 0 . 
ESTIMATION OF SOLUTION OF THE TASK 



(12) 
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Theorem 1. Let in (l)-(2) u 0 (x) < u + (0,x), xe R N ,g = xlT lp and 



1 



N 



P-m(p-\) p 



<0. 



Then in Q there is a global solution of (1), (2), for which we have the estimate u(t, x) < U + (t, x) , where 
u + (t,x) = u(t)f Q (^), (13) 
here u(t), f 0 (<^) given by formulas (3) and (11). 

Suggested that, /? > 1 , 0 < U 0 (x) < 1 , 0 < k(t)e C(0,°°) . For solution of this task is fair 

Theorem 2. Let 0< u 0 (x)< 1 xe R N . Then for the solution u(t, x) of the task (1), (2) in Q takes place an 



estimate 



u(t)(T Q +ty N e [4D,J <u(t,x)<e° (T 0 +ty"e 



jk(t)dt 



N „ [4Dt 



where U (?) -solution of the equation (1) without diffusion part: 



(14) 



U t — k(t)u(l-U^) ,i.e. u(t) 



and constant T n >l . 
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Study various properties of the solution of the Kolmogorov-Petrovsky-Piskunov-Fisher task with diffusion with 
reaction applied to the study of biological populations attracts the attention of many researchers, as it will be of interest 
from the point of view of mathematics as a nonlinear, which has an important application task. Therefore it's important 
generalizate this task on different cases. 

In particular, it has not been studied in case of a heterogeneous environment (diffusion coefficient is a function of 
the spatial variable), reaction coefficient depends on time or has a more complicated character. Consider the task of 
reaction diffusion of Kolmogorov - Fisher type 



du _ d 
dt dt 



du' 



dx 



du m 1 
dx 



+ k(t)u(l-un 



(16) 



u l (=0 = u Q (x),xe R N ,supw 0 (x) < +°o, 



(17) 



describing the process of diffusion-reaction in heterogeneous (anisotropic) environment with heterogeneity 
coefficient x" , where O^flG R . 

Propose an algorithm for constructing upper solutions of the Cauchy problem (16), (17) by the method of 
nonlinear splitting and etalon equations. 
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For the solutions of the Cauchy task (16), (17) fair following theorem: 
Theorem 3: Suppose that in (16) performed following conditions 

u(0,x) < Z(0,X) , where z(t,x) = u(t)-f(£) , 

u(t) is a solution of the equation U t — k(t)ll (1 — U ^ ) , 
P-1 

/(£) = a-b^- 1 



\x\ m(p-l) + n-l 

, where b — f" , a >£>, " — p 



^T{t)=\[m] m(p ~ X)+n ~ i dt. 

Then the problem (16), (17) is globally solvable and for the solution in the field of Q take place following 

estimation u(t,x) < z(t,x). 

I I" 

Presence of coefficient \x\ in (16) makes difficult of study the task (16), (17) in terms of building self-similar or 

some partial solutions, as well as numerical modeling. Turning to the new variables, assuming u(t,x) = w(t,(p(x)) 

and choosing <P(x) , reduce (16) to a form which in the case f5 > 1 previously more detail studied by many authors [1], 
[2], [5]. 

NUMERICAL SCHEMES AND METHODS OF SOLVING 

To solve numerically the task (1) - (2) in the £1 let' s build uniform grid W h by x with step h : 

co h = {x t = ih, h>0, i = 0,1,..., «, hn = b], 

and temporary grid 

m T = {tj = jT, t>o, y = o,i,...,m, tm = T}, t>o. 

Replace the task (18)-(20) by using balance method with implicit difference scheme and receive differential task 
with error 0(h 2 + T) . 

Now we present the results of numerical experiments. Iterative process is based on the method of Picard and 

Newton. 

Results of computational experiments show that all iterative methods are suitable for the constructed schemes. 
Below listed results of numerical experiments at k(t, x) = k(t) in a heterogeneous environment in two-dimensional case, 
for different values of parameters included in the equation (Fig.l): 
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re [0,1]; x 1 e [-6,6]; x 2 e [-6,6]; n Xi = 30 ; n Xj = 30 ; «, = 100 



Values of Parameters 



Results of Computational Experiment at 
the Initial Moment of the Time 



Results of Computational Experiment at 
the Final Moment of the Time 



m(p-V) + n -1 > 0 

m = 0.7, ,0 = 1.5 

n=3.1,p=4.1 

eps = 10" 3 

Average It 1=1, 43 
Average It2= 1 
Average It=2,43 





time(l) 



time (100) 



m(p-T) + n-l = 0 
m = 0.5, J3 = 1.5 

n=l,p=0.5 
eps = 10" 3 

Average It 1=2, 19 
Average It2= 1 ,5 1 
Average It=3,70 





time(l) 



time(lOO) 



m(p -l) + n-l< 0 
m = 0.1, £ = 1.5 

n=l,p=0.1 

eps = 10" 3 
Average Itl=2,02 
Average It2= 1 ,2 1 
Average It=3,23 





time(l) 



time(lOO) 



Figure 1: Results of Numerical Experiments in the Two-Dimensional Case 
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